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Abstract 

We use the Quasiparticle Random Phase Approximation (QRPA) and the Skyrme interactions 
SLy4 and SkM* to systematicahy calculate energies and transition strengths for the lowest 2"'" state 
in spherical even-even nuclei. The SkM* functional, applied to 178 spherical nuclei between Z = 10 
and 90, produces excitation energies that are on average 11% higher than experimental values, with 
residuals that fluctuate about the average by —35% + 55%. The predictions of SkM* and SLy4 have 
significant differences, in part because of differences in the calculated ground state deformations; 
SkM* performs better in both the average and dispersion of energies. Comparing the QRPA results 
with those of generator-coordinate-method (GCM) calculations, we find that the QRPA reproduces 
trends near closed shells better than the GCM, and overpredicts the energies less severely in general. 
We attribute part of the difference to a deficiency in the way the GCM is implemented. 

PACS numbers: 21.10.Pc, 21.60.Jz 
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I. INTRODUCTION 



Computer resources now allow nuclear Density Functional Theory (DFT), also called 
self-consistent mean field theory, to be applied systematically over the entire nuclear chart. 
Although it is usually applied to ground states, static DFT can be extended to treat exci- 
tations. A number of methods to do so have been developed, and it is not yet clear how 
accurate any of them is and how they compare to one another. In this paper we apply 
one of the methods, the Quasiparticle Random Phase Approximation (QRPA), to calculate 
the properties of the first excited 2'*' states in even-even nuclei. We assess strengths and 
weaknesses of the method for spherical nuclei and compare results with those of two other 
systematic studies that used different methods. We have previously carried out QRPA cal- 
culations for several isotope chains leading to the neutron drip line ^j] ; the same theory and 
codes are used here to cover essentially all known lowest-lying 2"*" states in spherical nuclei. 

Among the extensions of static DFT, the QRPA is one of the best justified theoreti- 
cally. It is the adiabatic limit of time- dependent DFT and as such may be derived from 
a variational principle. Time- dependent DFT preserves dynamical conservation laws, and 
in the adiabatic approximation is unambiguously defined by the density functional used for 
the static properties. Its major shortcoming as an approximation to many-body dynamics 
is that it works reliably only when ground state is nearly a single Slater determinant or a 
BCS-like condensate. 

Recently, two other systematic calculations of 2+ excitations that use very different ap- 
proximation schemes [2I, y] have appeared. Both of these studies constrain DFT by applying 
external quadrupole fields. This allows one to build a multi-configurational wave function or 
construct a collective potential-energy surface in deformation space. Reference 0], using the 
SLy4 energy-density functional [4] , treats many constrained configurations together to make 
a discrete-basis Hill- Wheeler approximation. Reference [sf, using the Gogny interaction [sl 
constructs a collective Hamiltonian in the 5-dimensional quadrupole-deformation space via 
the collective potential energy surface and a local treatment of the kinetic energy operator. 
We will refer to both as Generator- Coordinate Method (GCM) calculations. 

In addition to the approximation scheme, the quality of the energy-density functional 
also affects accuracy; most functionals are constructed to optimize their static properties. 
In this study we cannot hope to make better functional, but we can at least compare the 
quality of those that already exist. Here we will test two Skyrme functionals, SLy4 [3], and 
SkM* |6|, both in conjunction with "volume" pairing interactions. The results for SLy4 can 
be directly compared to those of Ref. [2], which uses the same functional^. We should note 
that a relativistic functional has already been applied in the QRPA to a subset of the nuclei 
we examine here 

II. IMPLEMENTATION OF QRPA 

Our QRPA and underlying Hartree-Fock-Bogoliubov (HFB) calculations assume spherical 
symmetry and employ a radial mesh with box boundary conditions. The HFB orbitals are 
calculated with a code based on the one described in Ref. |8|]. For the global survey here, we 



In our previous work we also considered the SkP functional [8j, but found that under certain conditions 
it lacks a true variational minimum . 
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give the box a radius of 16 fm and a radial mesh size of 0.05 fm. Although we can take the 
ordinary DFT functional from mass fits in the literature, the choice of a pairing functional 
is less clear cut. We have chosen to treat pairing through an ordinary contact interaction, 
— Vo5(ri — r2). After truncating the space of single-particle orbits (see below), we make 
an approximate fit of the pairing strength to experimental data. The quantities compared 
are the average HFB pairing gap and the 3-point difference A*^^-* of the experimental 
binding energies. The estimated strengths with the SLy4 functional are Vq = 200 MeV fm^ 
for neutrons and Vq = 240 MeV fm^ for protons ^. 

We use the matrix representation of the QRPA in a two-quasiparticle basis, as described 



in detail in Ref. [12|. We represent the HFB eigenf unctions in the canonical basis, and 
truncate them through limits on the occupation number (f ^) or, if the pairing gap is zero, 



on the single-quasiparticle energy 12]. We then construct the usual A and B matrices 
in the truncated two-quasiparticle basis. 

The method was applied to the Ca, Ni, and Sn isotopes for 0"^, 1~, and 2+ states in Ref. 
[l|. In this paper we apply the QRPA to many more nuclei than in Ref. [l|, but only consider 
the low-lying 2"*" states. Unlike the O"*" and 1~ excitations, the 2+ states are not affected by 
spurious modes in spherical nuclei; we can therefore afford to use a smaller space. In our 
prior work, for example, we used a box radius of at least 20 fm. Here we reduce the radius 
because the lowest 2+ states are well localized, even near the drip line. Because we do not 
have to be as careful about translational symmetry, we reduce the quasiparticle-energy cutoff 
in the HFB calculation from 200 MeV to 50 MeV. Finally, we increase the canonical-basis 
occupation- number cutoff — the smallest occupation number that canonical states included 
in the QRPA can have — from 10^^ to 10"'', and reduce the single-particle-energy cutoffs 
(used when pairing is absent) from 100 MeV to 30 MeV. The only way in which we must 
extend the basis is by increasing the cutoff in the single-particle angular momentum from 
21/2 (in our previous work) to 25/2. The reason is that we are dealing with heavier nuclei 
here. With all of these limits in place, our largest QRPA matrix, which we encountered in 
the very heavy nucleus ^^°Th, has dimension 6086. The space defined by all the cutoffs is 
large enough so that energies and B{E2)^ strengths do not change when we make it larger 
still. With the space size fixed, we then determine the strength Vq of the pairing interaction 
as described above. 



III. SELECTION OF NUCLEI AND THEIR QRPA EIGENMODES 

Our goal is to evaluate the performance of the QRPA in all the s phe rical nuclei with 
10 < Z < 90 for which experimental data have been tabulated in Ref. [ij]. That reference 
reports 507 nuclei with measured 2"*" excited states. We select spherical nuclei as follows: 
We use the HFB code ev8 jl5| to do constrained calculations on a grid of deformations, and 
drop the cases in which the lowest-energy configuration has nonzero deformation. For the 
remainder, we perform unconstrained calculations, starting from deformation values Q= 
and ±200 fm^. If any of these configurations evolves to a lower energy than the spherical 
one, we drop the nucleus. 



^ The finding that the extracted proton pairing strength is larger than the neutron value appears in other 
studies as well, e.g. Ref. The strengths for the SkM* functional are Vq — 170 MeV-fm'' for neutrons 



and Vq = 200 MeV-fm^ for protons. 
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FIG. 1: Particle-hole character of the lowest 2^ solutions. The histogram displays the quantity 
AA^ defined in Eq. (1) for 155 nuclei in the SLy4 data set (one of which we drop — see text). 
The values —2, 0, +2 correspond to excitations of hole-hole, particle-hole, and particle-particle 
character, respectively. 

A final requirement for selection is that the QRPA eigenvector correspond to a physical 
excitation. Thus, we exclude two nuclei in the SLy4 data because the eigenvalues are 
imaginary^. Other unphysical excitations are eigenvectors whose predominant components 
correspond to pair addition or removal, rather than excitation of the ground state. These 
components are present because the HFB quasiparticle vacuum does not have a well-defined 
number of particles. We exclude eigenvectors that have nearly pure two-particle-transfer 
character, using as a measure an approximate difference in particle number between the 
ground and excited states: 

^ = 2^:jXl^-Y^;){\-vl-v^. (1) 

Here is a canonical-basis condensate occupation probability, i,j label the quasiparticles 
in the two-quasiparticle basis, and X, F are the usual QRPA amplitudes. If |AA^| is large, 
we select another excitation or drop the nucleus from the table altogether. The final tables 
contain 155 spherical nuclei for the SLy4 functional and 178 for the SkM* functional. We 
note that the two tables have 129 out of a possible 155 nuclei in common. The imperfect 
overlap indicates that sphericity is not a very robust property of the DFT, unlike strong 
deformation. The data sets are posted with the calculated excitation properties in Ref. [l6j . 

Figured] shows a histogram of values of AA^ for the SLy4 data set. There is a well-defined 
peak at AA^ = with a spread of ±0.4. Beyond that, a plateau covers the range ±1.3. All 
these nuclei are included in the data set. Finally, there are 5 transitions that are clearly at 
the two-particle transfer limit AA^ ^ ±2. These nuclei are ^°'^®Ca, ^®Ni, ®'^Zr, and ^^^Sn. 
For all but one (^°Ca), an acceptable transition is present at an energy slightly higher than 
that of the lowest-energy solution and the nucleus is included in our survey. 



^ In principle, the stability of the spherical ground state demands that the eigenvalues all be real. However, 
the difference between the ways we truncate the HFB and QRPA calculations gives rise to occasional 
violations. 
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FIG. 2: Distribution of the summed {X + Y)'^ amplitudes, Eq. ([2]), for the first transitions in 
spherical nuclei, for the SLy4 functional. Values greater than 4 are consolidated in the last bin. 



We will also want to examine how well the QRPA handles "soft" nuclei, i.e. those which 
make large excursions from the sphericity. The size of the Y amplitudes is a good measure of 
softness. The enhancement of T-even transition rates due to the added ground state QRPA 
correlations is given roughly by the softness parameter 

C = T.i^^, + Y,,r. (2) 

This expression is designed to reduce to unity in the absence of the Y amplitudes and have 
the basic {X + Y^ structure characteristic of QRPA transition rates. A histogram for C 
appears in Figure [21 The softness factors range from 1 (no softness) to more than 4. Even 
the values in the middle of the histogram would make the Quasiparticle Tamm-DancofF 
Approximation invalid. The QRPA itself can have problems with very soft nuclei; if the 
Y amplitudes are comparable in size with the X amplitudes, then a phase transition to a 
ground-state of totally different character is nearby. 



IV. ENERGIES 

Figure [3] shows our predicted energies for the lowest 2+ states with the SLy4 functional; 
the predictions are plotted versus measured energies. The data span a range of more than 
an order of magnitude and the theory captures much of that variation. There are some cases 
where the errors are very large, however. Far off the line on the upper left hand side are 
nuclei with N = 40, Z = 34 — 40, shown by triangles. They are predicted to be spherical 
by the SLy4 functional but experimentally they appear to be deformed. On the lower side 
of the diagonal, the outlying cases are the isotopes of Sn (Z = 50) with neutron numbers 
N = 60,62,64 (inverted triangles). In this case, the SLy4 functional predicts the nuclei to 
be very soft with respect to quadrupole deformation. In fact, these 3 nuclei have the highest 
values of the softness parameter in the data set. 
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FIG. 3: Calculated energies for lowest 2^ states in spherical nuclei, plotted versus experimental en- 
ergies. The theoretical energies are from QRPA with the SLy4 energy functional. The experimental 
data are from Ref. 



14( 1 . For triangles, see text. 




-2 -1.5 -1 -0.5 0.5 1 1.5 2 



Re 

FIG. 4: Histogram of the quantity Re for the 155 nuclei in the SLy4 data set. The highest bin 
includes 3 nuclei having Re > 2. 

We can quantify the performance of calculations through the measure 

Re = ln(i?calc./-£^exp.) , (3) 

as is done in Ref. [2|. Figure H] shows a histogram of Re for the 155 nuclei in the SLy4 
data set. We see that theory tends to overpredict the excitation energy, but the error is 
systematic and the overall distribution is strongly peaked. To summarize performance in 



a single number, we can use the average Re for the data set. For spherical SLy4 nuclei, 
Re = 0.33, corresponding to a calculated energy about 40% higher than the experimental 
value. Another indicator of theory's performance is the width of the peak in the histogram. 
We define the dispersion of Re about the average as 

aE = ^J{Rl)-Rl. (4) 

A small value of aE implies that the theory tracks fluctuations in the data (though the 
average is not necessarily correct). For SLy4, the dispersion is aE = 0.51, corresponding 
to error bars of +66/ — 40 % in the energies. It should be noted that the peak is far from 
Gaussian, and the bounds icxs include more of the data set than would be the case for a 
Gaussian distribution. If we were to use the fraction of nuclei inside the bounds to define 
aE, its value would be 20% lower. 

It is interesting to see whether the accuracy of the energy prediction depends on other 
characteristics of the state. To examine this question, we use the value |AA^| to split the 
data set into two parts. Dividing the data set at |AA^| = 0.5 gives two roughly equal-size 
subsets. Table I shows the performance measures for the two subsets, and for reference 
the combined results discussed in the last paragraph. Performance in the high-|AA^| set is 
slightly poorer for both the average and the dispersion, but in our opinion, the differences 
are not large enough to warrant the use of |AA^| as a selection criterion beyond the most 
extreme cases. We also checked the dependence of accuracy on softness, dividing the data 
into subsets with C greater and less than 2. The results appear in lines 4 and 5 of the 
table. The separation strongly affects the averages, so that the high-softness set is actually 
more accurately described, but the dispersions are not significantly different. It is not totally 
surprising that states in soft nuclei are better described, as long as the nuclei are not too soft, 
because such states tend also to be more collective and RPA methods capture collectivity 
well. 

We next examine the performance of the SkM* functional. Here the protocol for selecting 
spherical nuclei gives a data set with 178 members. The set does not include the nuclei near 
®°Zr whose energies were poorly predicted by the SLy4 functional. Figure [5] shows the scatter 



TABLE I: Averages Re and standard deviations cte for measured 2"*" excitations with the SLy4 
and SkM* functionals and various cuts on the data set. The set labeled "common" consists of the 
nuclei that are spherical for both the SLy4 and the SkM* functionals. 



functional 


data set 


Number of nuclei 


Re 




SLy4 


all spherical 


155 


0.33 


0.51 




low |AiV| 


79 


0.29 


0.47 




high |AiV| 


78 


0.38 


0.54 




low softness 


106 


0.47 


0.48 




high softness 


49 


0.04 


0.44 




common 


129 


0.26 


0.40 


SkM* 


all spherical 


178 


0.11 


0.44 




low softness 


115 


0.27 


0.35 




high softness 


63 


-0.17 


0.45 




common 


129 


0.14 


0.38 
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plot of theoretical versus experimental excitation energies. There is more clustering along 
the diagonal than for SLy4, but outliers still exist on the lower side of the plot. These are 
all open-shell nuclei, and it is not as easy to characterize them as it was for SLy4. The worst 
cases are the nuclei ^^Ar, ^^Zn, and ^^^Th, shown as the inverted triangles in the figure. 
What they have in common is a very high QRPA softness, in the range C = 2 — 8. 

The R metrics for the SkM* data set are given on line 7 of the Table. The average is 
much better (lower both in the high- and low-softness set) than with SLy4, and there is also 
some improvement in the dispersion. Much of the improvement comes undoubtedly from 
the better selection of spherical nuclei, but part may be due to differences in the effective 
QRPA interaction that comes from the functional. One way to test that is to compare the 
performance on the 129 nuclei that the two data sets have in common. The results are labeled 
"common" in Table I. The averages still differ significantly — ^^(SLy4) — ^£;(SkM*) = 0.12 
— indicating that the functional plays a role beyond merely shaping the ground state. It 
remains to characterize the aspects of the functional that are responsible, but this is beyond 
the scope of our project here. 

The small- amplitude assumption of QRPA is best justified in magic and semimagic nuclei. 
Figures [6] and [7] show the very accurate QRPA results with SkM* in the vicinity of doubly 
magic nuclei in several isotope and isotone chains. Going from left to right in the figures, 
the magic numbers are {N,Z) = (28, 20), (28, 28), (40, 28), (50, 40), (82, 50), and (126,82). 
This last doubly magic nucleus is the one with the highest 2+ energy in all the chains. We 
obtain accurate results around doubly magic nuclei with SLy4 as well. 



V. TRANSITION STRENGTHS 



The 0+ 2+ transition strength for 88 of the 155 nuclei in the SLy4 spherical data set 
have been measured, and are compiled in Ref . |U] . Figure [8] shows a scatter plot comparing 
the SLy4 predictions of B{E2)^ with the experimental numbers. Four nuclei for which the 
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FIG. 5: Same as Figure O using the SkM* functional for the calculated energies. 
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FIG. 6: 2+ excitation energies in the vicinity of doubly magic nuclei. Displayed on either side of 
the doubly magic nucleus (A^, Z) are the adjacent even-even isotopes, (A^ — 2, Z) and {N + 2,Z). 
Theoretical energies (open circles) were calculated with the SkM* functional. Filled circles are 
experimental energies. 
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FIG. 7: Same as FigureOfor isotones (N, Z — 2), (N, Z + 2) adjacent to doubly magic nuclei {N, Z). 



theory seriously underpredicts the data are marked by triangles on the lower side of the 
scatter plot. The nuclei are '^^Ca, ^''Ar, ^^Ca, and ^^Sr, going from left to right on the plot. 
The nucleus ^^Sr is in the deformed {N, Z) ~ (40, 40) region; its large experimental transition 
strength confirms its deformed character. The failure of the theory for the Ca isotopes near 
= 20 can be explained by the absence in our calculation of deformed intruder orbits. 
Only one nucleus on the list, '^^Ar, does not have any obvious properties that would cause 
the QRPA to fail. 

On the upper side of the plot, we have marked the most prominent outlier, the nucleus 
^^°Po. Its measured strength is nearly an order of magnitude smaller than the calculated val- 
ues for both the SLy4 and the SkM* functionals. This discrepancy is one of the most puzzling 
in the survey, because deformation effects would only increase the transition strength. The 
predominant amplitude in the transition, with 93% of the normalization sum J2 X'^ — Y'^ = 1, 
has both quasiparticles in the lowest hg/2 states, corresponding to the shell- model transition 
{hg/2 — > (^9/2 ^9/2)''"^ • By itself, the pure shell-model transition has a strength 
of 454 e^fm^. The additional small X amplitudes'* in the QRPA increase the strength by 
factor of 3.4. An enhancement of the single-particle strength, often represented by effective 
shell-model charges, is nearly universal in quadrupole transitions. The nucleus ^^°Po is an 
exception, with a measured value 200 e^fm*, even smaller than the shell-model result with 
bare charges. 

To summarize the performance of the QRPA for transition strengths, we follow Ref. [2| 



The Y amplitudes play a minor role; the softness parameter C in Eq. ^ is close to unity (C — 1.09). 
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TABLE II: Averages Rq and standard deviations ctq of the distributions in Rq for the SLy4 and 
SkM* functionals. 









functional 






SLy4 


-0.32 


0.42 


SkM* 


-0.29 


0.53 



and define the residual of the transition matrix element logarithm Rq: 

(5) 



B{E2) Tcaic. /B{E2)] 



exp. 



Figure [9] shows a histogram of this quantity for the SLy4 functional. ^^''Po on the right 
clearly stands out as an anomalous case. The long tail on the left represents mainly deformed 
nuclei that are incorrectly predicted to be spherical. The average and standard deviations 
of Rq for both functionals are given in Table [TTl There is very little difference between the 
performance of the two functionals when it comes to transitions. 



VI. COMPARISON WITH GCM CALCULATIONS 

We refer to the GCM-based techniques of Refs. Q and Q as GCM-HiU- Wheeler (GCM- 
HW) and GCM-5-Dimensional-Collective-Hamiltonian (GCM-5DCH) respectively. The 
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FIG. 8: The same as Figure [3]but for B{E2)]. 
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FIG. 9: Histogram of Rq for the 155 nuclei in the SLy4 data set. The highest bin represents the 
nucleus ^^'^Po. 



GCM-HW study is in some ways the most like ours. It used the SLy4 functional, mak- 
ing the results directly comparable to our SLy4 results. And unlike the GCM-5DCH, it was 
applied to doubly magic nuclei. Its set of spherical nuclei is somewhat different from ours 
because of differing pairing interactions, so we use overlap of the two sets for comparison. 

The results are in Table III. The energy measures show some interesting differences. 
From the Re values we see that the GCM-HW gives energies about 40% higher than those 
of the QRPA. A possible explanation for the overprediction of energies by the GCM-HW 
in spherical nuclei appears in the Appendix A. The other interesting difference, which we 
cannot explain, is that the GCM-HW fluctuations track the experimental ones better than 
the QRPA. 

The B{E2)] metrics appear in the last two columns of the table. While the QRPA 
seriously underpredicts the transition strengths, the GCM-HW errs in the opposite direction. 
It is not surprising that the GCM wave function produces larger quadrupole matrix elements 
because it incorporates large deformation (unlike the QRPA), and because the J = and 
J = 2 states have a similar intrinsic structure. 

Unfortunately, because Ref. fsj used a different functional, it is more difficult to compare 
with our study; it is not possible to distinguish effects of the functional from those of the 
methodologies. We compare results nevertheless. As before, we use spherical nuclei analyzed 
in both studies to compute the /^-statistics; doubly magic nuclei are thereby excluded. We 
show in the first line of Table IIVI the /^-statistics of the QRPA with our better-performing 
functional, SkM*; in the second line are the results from Ref. 0] for the Gogny functional. 
We see that the QRPA (SkM*) has a somewhat better average energy. On the other hand, 
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TABLE III: Averages Re and standard deviations aE for 153 excitations produced with the 
SLy4 functional by our QRPA calculations and the GCM-HW calculations 

theory Re (Te Rq oq_ 

QRPA 0.33 0.51 -0.32 0.42 

GCM-HW 0.67 0.33 0.16 0.41 



TABLE IV: Same as Table IIIH with the SkM* functional for the QRPA and the Gogny functional 
for the GCM-5DCH. 



theory 


Re 


cte 


Rq 




QRPA (SkM*) 


0.10 


0.45 


-0.29 


0.51 


GCM-5DCH (Gogny) 


0.19 


0.43 


0.22 


0.27 



the fluctuation of the transition matrix elements is much better described by the GCM- 
5DCH (Gogny) treatment. 

We turn flnally to semimagic isotope chains. In Figure [TU] we show of our predictions 
(SkM*) and those of Ref. ^ (along with experimental data) for energies along 3 isotopic 
chains — Z = 28, 50, and 82 — and isotonic chains with the same magic numbers. Here our 
calculation is again systematically better than that of Ref. [2|, the results of which we omit 
from the plot for the sake of clarity. Our calculation and the GCM-5DCH calculation of 




FIG. 10: The lowest 2^ energies of even Ni, Sn, and Pb isotopes (panels a and b), and isotones 
with A'^ = 28, 50, and 82 (panel c). The open circles represent our results (SkM*), the squares the 
GCM-5DCII results of Ref. 0], and the filled circles experimental data. 
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Ref. [sl have complementary virtues in the isotopic chains. The nearly unchanging energies 
for Z = 50 (52 < iV < 80) and Z = 82 (112 < N < 124), used to motivate the idea of 
generalized seniority [l7|, have long challenged mean-field-based approaches. In Z = 50, 
for example, the GCM-5DCH reproduces this trend quite well up to = 70, but then 
incorrectly predicts a gradual increase in excitation energy as the closed shell at = 82 is 
approached. Our calculation, by contrast, produces lower energies than experiment around 
= 62 because the QRPA solutions are close to a transition to quadrupole deformation, 
but accurately reproduces the sharp jump at A^ = 82. In the Z = 82 chain we reproduce 
the energies well for 112 < N < 124, a bit better than does the GCM-HW, particularly near 
the closed shell. In the isotonic chains, by contrast, the results of the GCM-HW are pretty 
uniformly better than ours. Again, we cannot be sure how much of the difference is due to 
functionals and how much to methodology. Ref. [tI], which uses a relativistic functional to 
examine some of these same chains in the QRPA, obtains results that appear comparable 
to ours. 



VII. CONCLUSION 

We have used the QRPA to calculate the energies and E2 strengths of the lowest 2+ states 
in a wide range of even-even spherical nuclei, and compared the results with experiment for 
more than 150 energies and more than 80 strengths. We applied two functionals, SLy4 
and SkM*. On the whole, SkM* performed better than SLy4. For energies, our calculation 
with SkM* is comparable to the GCM-5DCH calculation of Ref. in spherical nuclei — 
better near closed shells, though perhaps not quite as good at midshell — and better than 
the GCM-HW calculation of Ref. Q. For 5(E2) t values, the GCM-5DCH calculations 
appear to be the best if doubly-magic nuclei are excluded. The QRPA is a small amplitude 
approximation, and its results are not disappointing, considering that limitation. 

We were surprised to find that the two functionals disagreed significantly on the question 
of which nuclei are spherical, an issue that arises at the mean-field level and has nothing to 
do with the QRPA. We plan to broaden our study soon to include deformed nuclei, allowing 
a more comprehensive comparison of approaches and functionals. 



VIII. ACKNOWLEDGMENT 

We thank A. Bulgac for discussions. This work was supported by the UNEDF SciDAC 
Collaboration under DOE grant DE-FC02-07ER41457. We used computers at the National 
Energy Research Scientific Computing Center at Lawrence Berkeley National Laboratory. 



APPENDIX A: THE RPA AND THE GCM 

In this appendix, we present a simple model that indicates why the GCM overpredicts 
the energy of the lowest 2+ excitation under conditions for which the RPA works well. The 
effect of the constraining field Q on the mean-field ground state |0) can be expressed in 
perturbations theory (that is, for small deformation) as 

i,)=io)-i:^^ (Ai) 
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where |0) is the spherical static-DFT state, \i) the RPA eigenstate, and Ei the corresponding 
RPA energy. In the GCH-HW approach, one first determines the polarized state \q). The 
state is then projected so as to have a well-defined angular momentum. The excitation 
energy is taken to be the expectation value of the Hamiltonian in the projected state. In 
Eq. (lAip . the projection ehminates |0) but keeps the relative amplitudes of the excited 
components in the above wave function. The expectation value of the Hamiltonian is then 
given by 



{E)q = ^ ^,|Q|op (A2) 

To get a quantitative estimate of error in the GCM-HW energy in this small-amplitude 
limit, we make a two-state approximation to the RPA spectrum, which in fact in spherical 
nuclei is often dominated by two strong transitions: an in-shell transition at a few MeV or 
less and the giant isovector quadrupole resonance at 10—15 MeV. The properties we need 
are their transition strengths and the ratio of their energies r = E1/E2. It is convenient to 
express the transition strengths as a fraction of the energy-weighted sum rule: 

_ E,{om' 

j:,e,{o\q\j)^- ^ 

With some simple algebra the energy expectation value {E)q can be written in terms of the 
lowest RPA energy Ei as 

The right hand side approaches the correct value {= Ei) as r — > 0. 

One case for which the small amplitude approximation (and therefore the RPA) should 
work well is the doubly magic nucleus ^°^Pb. In our QRPA calculation with the SLy4 
functional, we find an excitation energy of 5.0 MeV with a transition strength that amounts 
to 17% of the isoscalar sum rule. The giant resonance is located at E2 ~ 12 MeV, giving 
r 0.4. Inserting these numbers in Eq. (lA4p results in a GCM-HW energy that is 37% too 
large, and is close to the value obtained in Ref. 0, S2(GCM) = 6.8 MeV. This nucleus and 
others near closed shells, it should be noted, are not good ones for the GCM because of the 
size of r. Away from doubly magic nuclei r is smaller, resulting in a smaller GCM error for 
a given /i. 

It might be possible to improve the GCM by adding an additional quadrupole field that 
would distinguish the giant resonance from the low-lying collective states. 



APPENDIX B: SUM-RULE FRACTION 

It has been known for a long time that the 2"'" excitation energies and the associated 
B{E2)'\ values have a strong inverse correlation [l8|. and it is interesting to see how well 
the theoretical results reproduce this behavior. The product E2B{E2)] can be expressed as 
a fraction of an energy-weighted sum rule, giving it the A-dependence of the sum rule. A 
form that is often used is Z^y4~", where a = 1/3 corresponds to the isoscalar sum rule for 
a liquid drop^^. We show in Table IVl the average of the scaled product 

Z^A-^I^ ^^^^ 
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TABLE V: Average sum rule fraction (Eq. ()Bip ) and variance 



source 


data set 


Number of nuclei 


S = exp{Rs) 


0-5 


experiment 


all 


328 


5.4 10"^ 


0.46 


experiment 


spherical 


99 


4.7 10"^ 


0.63 


DFT 


spherical 


99 


2.8 10-^ 


0.95 



obtained from averaging Rs = log(S'). The variance of Rs is also given in the table. The first 
line shows the values obtained from the experimental data [14]. The second line contains 
only the nuclei predicted to be spherical with the SkM* functional. One sees that the average 
over the spherical nuclei is close to the global average, indicating that the sum-rule fraction 
is insensitive to deformation. The variance is larger among the spherical nuclei, however. 
The third line shows the results obtained in the QRPA with the SkM* functional. The 
theoretical average value is significantly smaller than the empirical value. Since the fraction 
of the sum rule is a measure of the collectivity, it may be that our functional does not have 
sufficient collectivity. The sum-rule fraction might be sensitive to the strength of the pairing 
interaction or some other aspect of the functional, and it would be interesting to investigate 
that issue further. 

The table also shows that the fluctuation in the sum-rule fraction, as, is significantly 
larger in the DFT calculation than in experiment. Shell structure may need to be suppressed 
in some way, perhaps by increasing the pairing strength. 
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